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Abstract 

Dijet angular distributions of photoproduction events in which a D*^ meson is 
produced in association with one of two energetic jets have been measured with 
the ZEUS detector at HERA, using an integrated luminosity of 120 pb _1 . Dif- 
ferential cross sections as a function of the angle between the charm-jet and the 
proton-beam direction in the dijet rest frame have been measured for samples 
enriched in direct or resolved photon events. The results are compared with pre- 
dictions from leading-order parton-shower Monte Carlo models and with next-to- 
leading-order QCD calculations. The angular distributions show clear evidence 
for the existence of charm originating from the photon. 
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1 Introduction 



High-energy collisions at the HERA ep collider between a quasi-real photon and a proton 
provide an effective source of photoproduction processes. Jets with high transverse energy 
and/or charm (c) quarks produced in such processes can be described within quantum 
chromodynamics (QCD) in two sub-classes: direct processes are those in which the photon 
couples as a point-like object in the hard scattering; resolved processes are those in which 
the photon acts as a source of incoming partons, one of which participates in the hard 
interaction. Both processes can lead to two jets in the final state. Samples enriched in 
direct and resolved photon events can be identified using the variable x° bs P, which is 
the fraction of the photon's momentum contributing to the production of the two jets. 

Inclusive cross sections for photoproduction of D* ± (2010) mesons as well as cross sections 
for "charm dijet" events, in which the D* is observed in events with two energetic jets, 
have been previously reported |2j. Differential cross sections of the D* and associated 
dijet system are larger than next-to-leading-order (NLO) QCD predictions [3] at low x° bs , 
but are in agreement at high x° bs . The data were also compared to predictions of leading- 
logarithmic parton-shower Monte Carlo (MC) models. According to these comparisons, 
about 60% of the events can be attributed to the direct photon-gluon-fusion (PGF) process 
7# — > cc, illustrated in Fig. la). The MC models predict that most of the resolved photon 
events come from charm excitation of the photon (Figs, lc) and Id)) rather than from 
the gg — > cc process (Fig. lb)). The aim of this analysis is to determine the dominant 
mechanisms for charm dijet photoproduction in both direct and resolved photon processes. 

Measuring the angular distribution of the outgoing jets allows the dominant subprocesses 
to be determined and the MC predictions to be tested, as was done previously j3] for 
inclusive dijet events. This study showed that the differential cross-section da/d\ cos#*|, 
where 9* is the angle between the jet-jet axis and the proton beam direction in the dijet rest 
frame, is sensitive to the spin of the propagator in the hard subprocess. In direct photon 
processes, in which the propagator in the leading-order (LO) QCD diagrams is a quark, 
the differential cross section rises slowly towards high | cos 6 1 * | values (da/d\cos8*\ ~ 
(1 — | cos6*\)~ 1 ). In resolved photon processes, the gluon propagator is allowed at LO and 
dominates over the quark propagator due to the stronger gluon-gluon coupling compared 
to the quark-gluon coupling. In this case the cross section rises steeply when | cos 6* \ 
increases (da/d\ cos6>*| « (1 — | cos9*\)~ 2 ). Similar results have been reported in photon- 
photon collisions [Sj. 

If most of the resolved-photon charm dijet events are produced as a result of charm from 
the photon, a gluon-exchange contribution, as seen in Fig. Id), should dominate. This 
results in a steep rise of the cross section towards high | cos 6* | values. The other diagrams 
of Fig. 1 involve quark exchange and thus should not show such a sharp rise. If one of 
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the jets is explicitly tagged as a charm jet, the sign of cos 6* can be defined. If the charm 
originates from the photon, the charm jet generally lies in the photon hemisphere. 



2 Experimental conditions 

The analysis was performed using data collected with the ZEUS detector at HERA 
during 1996 - 2000. In this period, HERA collided electrons or positrons with energy 
E e = 27.5 GeV and protons with energy E v = 820 GeV (1996 - 1997) or E v = 920 GeV 
(1998 - 2000), corresponding to integrated luminosities of 38.6 ± 0.6 and 81.9 ± 1.8 pb -1 
and to centre-of-mass energies -y/i = 300 GeV and s/s = 318 GeV, respectively. This 
data sample is about a factor of three larger than that used for the previous charm dijet 
analysis [2]. 

A detailed description of the ZEUS detector can be found elsewhere A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [7j, which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD 
consists of 72 cylindrical drift chamber layers, organized in 9 superlayers covering the 
polar-angle 1 region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks is a{p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T (p T in GeV). 

The high-resolution uranium-scintillator calorimeter (CAL) jH] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a(E)/E = 0.18/ \/E for electrons 
and a(E)/E = 0.35/ \[E for hadrons (E in GeV). 

The luminosity was measured from the rate of the bremsstrahlung process e + p — > e + 7p, 
where the photon was measured in a lead-scintillator calorimeter [Hj placed in the HERA 
tunnel at Z = — 107 m. 



The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as rj = — In (tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. 
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3 Event selection 



Photoproduction events were selected with a three-level trigger [HIE]]. The inclusive 
photoproduction sample was defined by requiring a reconstructed vertex and no scattered 
electron or positron found in the CAL, thus restricting the photon virtuality, Q 2 , to be 
below 1 GeV 2 , with median Q 2 m 3- 10 -4 GeV 2 . The photon-proton centre-of-mass energy, 
W, was restricted to the range 130 < W < 280 GeV. The latter was measured using the 
Jacquet-Blondel [TTj estimator Wm = ^ Ay^E e E p ^ where t/jb = ^(Ei — pz,i) /2E e , the 
sum runs over all CAL cells and pz,i is the Z component of the momentum vector assigned 
to each cell of energy Ei. Jets were reconstructed with the cluster algorithm jT^j in its 
longitudinally invariant inclusive mode fH] . The events were required to have at least two 
jets 2 with pseudorapidity \r]^ ct \ < 2.4 and transverse energy Ej? > 5 GeV. The measured 
jet energies as well as Wjb were corrected for energy losses in inactive material in front 
of the CAL, using the MC simulation. 

The D* mesons were reconstructed using the mass-difference technique applied to the 
decay chain 3 D*^ 1 — ► D°n^ — > K^^ir^. Tracks in the CTD with opposite charges and 
transverse momenta pt > 0.5 GeV were combined in pairs to form D° candidates. Kaon 
and pion masses were assumed in turn for each track to calculate the pair invariant mass, 
M(Ktt). A third track, ns, assumed to be the "soft pion" from the D* decay, with 
Pt > 0.15 GeV and a charge opposite to the kaon, was added to form a D* candidate. 
Events with a mass difference AM = M{Kim s ) - M(Kir) in the range 0.1435 < AM < 
0.1475 GeV around the nominal value [Hj and the range 1.81 < M(Kn) < 1.92 GeV 
around the D° mass were called D* candidates. To suppress combinatorial background, 
a cut p®* / E T >10 ° > 0.15 was applied j2], where E T >10 ° is the transverse energy measured 
in the CAL outside a cone of 9 = 10° in the forward direction. The reconstructed D* 
mesons were required to have p® > 3 GeV and pseudorapidity in the range |?7 D *| < 1.5. 

These cuts ensure that the events lie in a well understood acceptance region of the detector. 



4 Jet kinematic variables 



Samples enriched in direct and resolved photon events were separated by a selection on 
the variable 



X 7 



S jets (4 ct e"^ 

2^Ee 



2 The fraction of events with more than two jets is 11%. 

3 Throughout this article, D° refers to both D° and D°. 
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where yE e is the initial photon energy and the sum is over the two jets with the highest 
Ej, . The selection of x° bs > 0.75 and x° bs < 0.75 yields samples enriched in direct and 
resolved photon processes, respectively. 

A complementary variable is 



S jets (4 Gt e" JCt ) 



which is the fraction of the proton's momentum contributing to the production of the two 
jets. 

The dijet scattering angle, 9*, is reconstructed using 



/ jetl _ „jet2\ 

cos0* = tanhf^ ^— J. (1) 

In the simple case in which two jets are back-to-back in the transverse plane and have 
equal transverse energies, the dijet invariant mass is given by Mjj = 2E i ^ t / ^/l — | cos6**| 2 . 
Therefore, for a given My, events with high values of | cos6**| have lower E^ . In order to 
study the | cos^*! distribution up to | cos9*\ = 0.83 without bias from the E 1 ^ 1 cut, Mjj 
was required to be above 18 GeV. 

A cut on the average longitudinal boost, fj = (^J ctl + r]i ct2 )/2, of |^| < 0.7 was applied. 
This selection limits rf ct to |77j ct | < 1.9 and removes the bias caused by the explicit cuts 
on rf et It also reduces the bias caused by the cut on \rj D *\ < 1.5 while retaining a 
sufficiently large number of events. Monte Carlo studies show that the residual distortion 
due to the \r] D *\ cut is small and confined to the extreme bins of the cos#* distribution. 

These cuts ensure that any features seen in the measured distributions can be attributed 
to the dynamics of the hard scattering processes. 



5 Models and QCD calculations 

The MC simulation programs PYTHIA 6.156 jUj and HERWIG 6.301 [H] were used to 
model the final states. The PYTHIA and HERWIG simulations use on-shell LO matrix 
elements for charm photoproduction processes. Higher-order QCD effects are simulated 
in the leading-logarithmic approximation with initial- and final-state radiation obeying 
DGLAP evolution Coherence effects from soft-gluon interference are included. The 
parton density functions (PDF) CTEQ5L [HI for the proton and GRV-G LO [Uj for the 
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photon were used. The LO direct and resolved photon processes were generated propor- 
tionally to their predicted MC cross sections, using charm- and beauty-quark masses of 
m c = 1.5 GeV and = 4.75 GeV, respectively. Fragmentation into hadrons is simu- 
lated in HERWIG with a cluster algorithm [2D] and in PYTHIA with the Lund string 
model |2U- 

Samples of MC events larger than the dataset were produced. To calculate the acceptances 
and to estimate hadronisation effects, the events were passed through the GEANT 3.13- 
based [22] simulation of the ZEUS detector and trigger. They were reconstructed and 
analysed by the same program chain as the data. Samples corresponding to different 
data taking conditions were generated in proportion to their luminosities. For PYTHIA, 
in addition to the D* decay chain used for this analysis, background events that arise 
from other D*^ decay modes or similar decay modes of other charm mesons were also 
simulated. 

The MC event generator CASCADE 1.00/09 j2H! simulates heavy-quark photoproduction 
in the framework of the semi-hard or /ct-factorisation approach [21] . The matrix element 
used in CASCADE is the off-shell LO PGF process. The CASCADE initial-state radiation 
is based on CCFM evolution [23], which includes in the perturbative expansion the ln(l/x) 
terms in addition to the InQ 2 terms used in DGLAP evolution. To simulate final-state 
radiation, CASCADE uses PYTHIA 6.1 and the fragmentation into hadrons is simulated 
with the Lund string model. The cross section is calculated by convoluting the off-shell 
PGF matrix element with the unintegrated gluon density of the proton obtained from the 
CCFM fit to the HERA F2 data, by fixing most of the free parameters j2H|. Although 
the CASCADE matrix element corresponds to the off-shell PGF direct photon process 
only (Fig. la)), resolved photon processes are reproduced by the CCFM initial-state 
radiation [2B]. 

The NLO QCD calculations of differential cross sections for photoproduction of charm di- 
jet events in the HERA kinematic region are available [3] in the fixed-order (FO) scheme. 
The PDF parameterisations used were CTEQ5M1 ^H] for the proton and AFGHO [21] for 
the photon. The factorisation scales of the photon and proton PDFs, fip, and the renor- 
malisation scale, /j,r, used for the calculation were set to fip — I^r — m T = a/ rn? c + (pj) , 
where (p^) was set to the average p\ of the charm quark and antiquark. The charm frag- 
mentation into D* was performed using the Peterson fragmentation function [2H] with an 
e parameter of 0.035 [29] . 

In all cases, the fraction of c quarks fragmenting into a D* was assumed to be 0.235 [HD] 
and a charm quark mass of m c = 1.5 GeV was used. 
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6 Results 



The AM distribution for dijet events in the D° signal region shows a clear D* signal. 
The analysis is based on 1092 ± 43 D*± mesons found in the 0.1435 < AM < 0.1475 GeV 
region over a background of 328 events. The signal has similar characteristics as that in the 
previous ZEUS publication [2] except that the signal to background ratio has improved by 
a factor of three due to the tighter cuts (see Sections 3 and 4) used here. The background 
was determined from the AM distribution for wrong-charge combinations, where the 
tracks forming D° candidates had the same charge and the tts had the opposite charge. 

The number of events in each bin of the measured variables was extracted by performing 
a bin-by-bin wrong-charge background subtraction. To obtain differential cross sections, 
each value was then multiplied by a correction factor proportional to the ratio of generated 
to reconstructed events from the PYTHIA MC simulation. The measured cross sections 
are the luminosity-weighted average of the cross sections at the centre-of-mass energies 
y/i = 300 GeV and y'S = 318 GeV. 

The systematic uncertainties were determined by adding the contributions from several 
sources in quadrature. The largest contributions were associated with the cuts on W and 
with the difference between the correction factors evaluated using HERWIG rather than 
PYTHIA. The uncertainties due to the knowledge of the CAL energy scale (±3%) are 
highly correlated between bins and are therefore shown separately. Statistical uncertain- 
ties dominate over systematic ones in most bins. The measured cross sections and their 
uncertainties are given in Tables 1-4. 

The differential cross section as a function of £° bs is shown in Fig. 2. The peak at high 
values of x° bs indicates a large contribution from direct photon processes, but there is 
also a sizeable contribution from resolved photon processes at lower x° bs values. Figure 3 
shows the differential cross section as a function of x° bs . The x° bs range of the data 
is concentrated in the region 0.0055 < x° bs < 0.044, where the proton PDFs are well 
determined. 

Figure 4 shows the differential cross sections as a function of | cos 6* | separately for the 
resolved-enriched (x° bs < 0.75) and direct-enriched (x° bs > 0.75) samples. The cross 
section for the sample enriched in resolved photons exhibits a more rapid rise towards high 
values of | cos 8* | than does the cross section for the sample enriched in direct photons. 
Consequently, the LO subprocess gg — > cc (Fig. lb)), with g-exchange in the t channel, 
cannot be the dominant resolved photon process for charm dijet events. This observation 
suggests a large gluon-exchange contribution originating from a charm-excitation process. 

The | cos#*| distributions of Fig. 4 are similar in shape to the previously reported dijet 
angular distributions [3] , which did not require the presence of charm. In those analyses, 
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only the absolute value of cos 6* was determined. In the present study, the two jets were 
distinguished by associating the D* meson to the closest jet in rj — (f) space. The associated 
jet is defined to be the jet with the smallest R { = yjjr^ - r] D *) 2 + - (p D *) 2 ; 
(i = 1,2) and with R < 1, where (ft ct (4> D *) is the azimuthal angle of the jet (D*) in the 
laboratory frame. Calling this ll D* jet" jet 1 in Eq. (1), the rise of da/dcosO* can be 
studied separately for the photon and proton directions. Figure 5 shows the differential 
cross sections as a function of cos 9* for the resolved- and direct-enriched samples. Events 
that did not satisfy the requirement R < 1 for at least one of the two jets (8.7% for 
x° bs < 0.75 and 1.1% for x° bs > 0.75) were not included in these cos^* distributions. The 
PYTHIA estimation of the contribution of the direct process to the resolved-enriched 
sample, x° bs < 0.75, and the resolved process to the direct-enriched sample, x° bs > 0.75, 
are also indicated. 

Direct photon events originating from the dominant g-exchange process jg — > cc (Fig. la)) 
should have a distribution symmetric in cos9*. The angular distribution of direct- 
enriched events (x° bs > 0.75) exhibits a slight asymmetry, which can be explained by the 
feedthrough from resolved photon processes near cos 9* = —1, as predicted by PYTHIA 
(Fig. 5b)). 

The sample enriched in resolved photons (Fig. 5a, c)) exhibits a mild rise in the proton 
hemisphere towards cos9* = 1, consistent with expectations from quark exchange. In 
contrast, they have a strong rise towards cos#* = —1, i.e. in the photon direction, 
consistent with a dominant contribution from gluon exchange. For the latter case, the 
charm quark emerges in the photon hemisphere (Fig. Id)). Gluon-exchange diagrams 
with this topology can only come, at LO, from the processes c 1 'g p — > eg and c 7 g p — > cq, 
where the superscripts refer to an origin in either the photon or proton. The partonic 
cross sections for these 2^2 subprocesses are highly asymmetric in cos 9* and show 
a steep rise towards the photon direction, while the subprocess gg — > cc (Fig. lb)) is 
symmetric in cos^*. This observation suggests that the source of the LO gluon-exchange 
contribution as seen in Figs. 4a) and c) is charm originating from the photon. This is 
consistent with the MC prediction |2] that most of the resolved photon contribution to 
charm dijet events at HERA is due to charm originating from the photon. 

7 Comparisons with theoretical predictions 
7.1 Comparison with MC predictions 

Figures 2-5 compare the distributions of the data with those of the MC simulations 
PYTHIA, HERWIG and CASCADE. For PYTHIA and HERWIG, the predictions are 
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normalised to the data with normalisation factors shown in brackets within the figures. 
For a shape comparison, the prediction for CASCADE is shown in Fig. 2 normalised to the 
data. Since there is a hope that higher-order corrections to fcj-factorised calculations 
might be smaller than those to LO parton-shower calculations using DGLAP evolution, 
the absolute predictions from CASCADE for the differential cross sections are shown in 
Figs. 3-5. 

The shapes of all data distributions are well reproduced by PYTHIA. The HERWIG 
predictions give an adequate description of the shapes in the data, although the rise in 
the cross section as a function of cos 9* at low x° bs is stronger in the data, particularly in 
the photon direction. There is a sizeable contribution from a resolved photon component 
in both PYTHIA (35%) and HERWIG (22%). Fitting the MC distributions to the data, 
allowing the resolved and direct photon contributions to vary independently, results in a 
resolved contribution of 46% for PYTHIA and 30% for HERWIG. The fraction of charm 
dijet events that originates from beauty production is predicted to be ~ 10% by PYTHIA 
and ~ 6% by HERWIG. The shape of the beauty component is similar to that of the 
overall distributions. 

The x° bs distribution of CASCADE, normalised to the data, gives a larger contribution 
at high x° bs and a smaller contribution at low x° bs (Fig. 2a)). The absolute cross section 
predictions for CASCADE, shown in Figs. 3-5, are larger than the data by around 30%. 
This difference is concentrated in the region x° bs > 0.75 and cannot be accounted for by 
a variation of m c : changing m c to 1.3 and 1.7 GeV gave a deviation in the prediction of 
±10%. However, the CASCADE prediction reproduces the shape in x° bs . The angular 
distributions are well described for x° bs < 0.75, although CASCADE underestimates 
the data in the proton direction (Fig. 5c). For x° bs > 0.75 (Fig. 5d), the prediction 
overestimates the data in all regions of cos 6*, although the shape is described reasonably 
well. 

7.2 Comparison with NLO QCD predictions 

The differential cross sections of Figs. 2-5 have been compared to the NLO FO calcu- 
lation Pj. The uncertainties in the NLO calculation, shown as the shaded area, come 
from the simultaneous variation of m c between 1.3 and 1.7 GeV and /i# between m-r/2 
and 2rriT- Changing the photon PDF parameterisation from AFGHO to GRVHO |3*2lll9j . 
as well as varying [ip of the photon and proton PDFs between m^/2 and 2mr, produce 
small effects (< 5%) on the NLO predictions. 

The differential cross sections predicted by the FO calculation were corrected for hadroni- 
sation effects. For each bin, the partonic cross section was multiplied by a hadronisation 
correction factor, Chad = <c r ° n V^MC ° nS ; which is the ratio of the MC 

cross sections after 
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and before the hadronisation process. The value of Chad wa s taken as the mean of the 
ratios obtained using HERWIG and PYTHIA. Half the spread between the two MCs was 
added in quadrature to the uncertainty in the NLO calculation. The deviation of C*h a d 
from unity is typically below 20% (see Tables 1-4). 

Figure 2b) shows a comparison for the differential cross section in x° bs . To minimise 
the large migration effects at x° bs > 0.75 due to hadronisation, a wider bin than that of 
Fig. 2a) was used. Migrations to low x° bs are small. The cross section can have a low 
x° bs contribution at NLO due to three-parton final states in which one of the partons is 
treated as a photon remnant. However, the low x° bs tail of the NLO cross section is below 
the data 0. For x° bs > 0.75, the data are well described by the NLO prediction. 

The differential cross section as a function of a;° bs is compared in Fig. 3b) with the NLO FO 
calculation. The NLO prediction is in reasonable agreement with the data. As expected 
from the x° bs comparison, the NLO prediction for the resolved-enriched x° bs distribution 
(not shown) is too low, but the shape is well reproduced. 

Figures 4c-d) and 5c-d) compare the charm dijet angular distributions to the NLO cal- 
culation. For high x° bs (Figs. 4d) and 5d)), the NLO prediction gives a good description 
of the data. For low x° bs (Fig. 4c)), the NLO prediction is significantly below the data. 
In Fig. 5c), the NLO predicts a lower cross section than the data in both proton and 
photon directions. The shapes of the | cos 9* | and cos 9* distributions are reasonably well 
described by the NLO predictions. 

8 Conclusions 

The differential cross sections as a function of cos 9* for charm dijet photoproduction 
events ( median Q 2 m 3- 10~ 4 GeV 2 ) have been measured in the kinematic range 
130 < W < 280 GeV, Q 2 < 1 GeV 2 , V t* > 3 GeV, \r] D *\ < 1.5, > 5 GeV and 
|^jet| < 2.4. The cuts on the dijet invariant mass, Mjj > 18 GeV, and on the average jet 
pseudorapidity, \fj\ < 0.7, select an Mjj and \fj\ region where the biases from other kine- 
matic cuts are minimised. The distributions have been measured separately for samples of 
events enriched in resolved (x° bs < 0.75) and direct (x° bs > 0.75) photon processes. The 
angular dependence for the two samples is significantly different, reflecting the different 
spins of the quark and gluon propagators. The cross section rises faster with increasing 
| cos 9* | for resolved photoproduction, where processes involving spin-1 gluon exchange 
dominate, than for direct photoproduction, where processes involving spin-1/2 quark ex- 
change dominate. 

The shapes of the measured differential cross sections are well reproduced by PYTHIA. 
Except for the angular distributions at low x° bs , HERWIG gives an adequate description of 
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these shapes. The predictions of CASCADE describe the data at low x° bs in both shape 
and normalisation. For high x° bs , the prediction significantly overestimates the data, 
but gives a reasonable description of the shapes. The shapes of the measured angular 
distributions are approximately reproduced by the NLO FO predictions. The absolute 
cross sections predicted by the NLO FO calculation reproduce the data for the sample 
enriched in direct photons but are below the data for the sample enriched in resolved 
photons. 

Associating the D* meson with one of the jets allows the sign of cos 6* to be defined. In 
all cases, the cos8* distributions show a mild rise towards | cos6*\ = 1, as expected from 
quark exchange, except for the resolved-enriched sample in which the cross section rises 
steeply in the photon direction (cos 8* = —1), as expected from gluon exchange. This 
observation indicates that most of the resolved photon contribution in LO QCD charm 
production is due to charm originating from the photon, rather than to the competing 
resolved photon process gg — > cc. This demonstrates that charm originating from the 
photon is the dominant component in the resolved photoproduction of dijet events with 
charm. 
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~,obs 
X 7 


bin 


dcr/dx° hs 


A s tat 


A S y S t 


A ES 


(nb) 


Chad 


0.250, 


0.375 


0.115 


±0.029 


+0.037 
-0.017 


+0.004 
-0.001 




0.941 ± 0.040 


0.375, 


0.500 


0.196 


±0.055 


+0.064 
-0.034 


+0.022 
-0.005 




0.950 ± 0.004 


0.500, 


0.625 


0.407 


±0.082 


+0.101 
-0.086 


+0.032 
-0.029 




1.006 ± 0.008 


0.625, 


0.750 


1.011 


±0.102 


+0.073 
-0.169 


+0.093 
-0.112 




1.285 ± 0.050 


0.750, 


0.875 


2.000 


±0.147 


+0.159 
-0.105 


+0.150 
-0.159 






0.875, 


1.000 


1.727 


±0.122 


+0.243 
-0.080 


+0.052 
-0.089 






0.750, 


1.000 


1.864 


±0.096 


+0.145 
-0.066 


+0.101 
-0.124 




0.851 ± 0.041 



Table 1: Measured cross sections as a function ofx° . The statistical, systematic 
and jet energy scale (Aes) uncertainties are shown separately. The multiplicative 
hadronisation correction applied to the NLO prediction is shown in the last column. 
The uncertainty shown for the hadronisation correction is half the spread between 
the values obtained using the HERWIG and PYTHIA models. 



x[ 



obs 



bin 



da/dx° p 



obs 



A s 



A SyS t 



A 



ES 



(nb) 



had 



0.0055, 
0.0110, 
0.0165, 
0.0220, 
0.0275, 
0.0330, 
0.0385, 



0.0110 
0.0165 
0.0220 
0.0275 
0.0330 
0.0385 
0.0440 



23.28 

36.90 

30.72 

18.55 

7.84 

3.21 

2.37 



±2.47 
±2.96 
±2.57 
±2.05 
±1.50 
±0.77 
±0.69 



+2.70 
-2.58 

+1.95 
-2.74 

+2.91 
-2.08 

+0.97 
-1.18 

+0.33 
-1.88 

+0.56 
-0.39 

+0.55 
-0.31 



+0.62 
-0.95 

+2.53 
-3.05 

+1.90 
-2.20 

+1.61 
-1.49 

+0.58 
-0.37 

+0.02 
-0.41 

+0.41 
-0.02 



0.799 ± 0.040 
0.910 ± 0.031 
0.953 ± 0.027 
0.985 ± 0.017 
0.984 ± 0.019 
1.020 ± 0.047 
1.022 ± 0.012 



Table 2: Measured cross sections as a function of x° . For further details, see 
the caption to Tabled 
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cos 6 


*| bin 


dajd\ cos#* 


A s tat 


A S y S t 


A ES (nb) 


Chad 


x° bs < 0.75 


0.00000, 


0.10375 


0.056 


±0.034 


+0.022 
-0.022 


+0.014 
-0.005 


1.007 ± 0.014 


0.10375, 


0.20750 


0.040 


±0.027 


+0.028 
-0.010 


+0.011 
-0.003 


1.099 ± 0.003 


0.20750, 


0.31125 


0.126 


±0.041 


+0.022 
-0.026 


+0.011 
-0.014 


1.072 ± 0.026 


0.31125, 


0.41500 


0.114 


±0.032 


+0.032 
-0.025 


+0.015 
-0.005 


1.099 ± 0.048 


0.41500, 


0.51875 


0.280 


±0.062 


+0.055 
-0.051 


+0.027 
-0.021 


1.107 ± 0.041 


0.51875, 


0.62250 


0.300 


±0.069 


+0.095 
-0.059 


+0.009 
-0.050 


1.101 ± 0.029 


0.62250, 


0.72625 


0.536 


±0.088 


+0.031 
-0.138 


+0.005 
-0.045 


1.145 ± 0.014 


0.72625, 


0.83000 


0.732 


±0.108 


+0.087 
-0.155 


+0.053 
-0.036 


1.115 ± 0.018 


if > 0.75 


0.00000, 


0.10375 


0.277 


±0.055 


+0.049 
-0.038 


+0.012 
-0.030 


0.923 ± 0.069 


0.10375, 


0.20750 


0.401 


±0.065 


+0.037 
-0.064 


+0.044 
-0.030 


0.919 ± 0.044 


0.20750, 


0.31125 


0.471 


±0.063 


+0.045 
-0.080 


+0.020 
-0.039 


0.910 ± 0.052 


0.31125, 


0.41500 


0.390 


±0.070 


+0.055 
-0.036 


+0.023 
-0.028 


0.906 ± 0.068 


0.41500, 


0.51875 


0.584 


±0.082 


+0.048 
-0.047 


+0.024 
-0.034 


0.876 ± 0.056 


0.51875, 


0.62250 


0.636 


±0.089 


+0.044 
-0.040 


+0.017 
-0.030 


0.863 ± 0.061 


0.62250, 


0.72625 


0.810 


±0.098 


+0.094 
-0.026 


+0.059 
-0.041 


0.832 ± 0.036 


0.72625, 


0.83000 


0.922 


±0.126 


+0.105 
-0.090 


+0.026 
-0.046 


0.756 ± 0.013 



Table 3: Measured cross sections as a function of |cos#*| for x° < 0.75 and 
x° bs > 0.75. For further details, see the caption to Tabled 
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cos#* bin 


da Jd cos 9* 


Astat 




A ES (nb) 


Chad 






™obs 

"S 


< 0.75 






-0.830, -0.664 


0.471 


±0.072 


+0.065 
-0.077 


+0.034 
-0.023 


1.063 ± 0.008 


-0.664, -0.498 


0.198 


±0.036 


+0.043 
-0.025 


+0.006 
-0.018 


1.065 ± 0.023 


-0.498, -0.332 


0.111 


±0.028 


+0.028 
-0.014 


+0.012 
-0.007 


1.084 ± 0.029 


-0.332, 0.000 


0.032 


±0.011 


+0.009 
-0.010 


+0.006 
-0.003 


1.056 ± 0.0004 


0.000, 0.332 


0.043 


±0.015 


+0.009 
-0.007 


+0.006 
-0.004 


1.105 ± 0.061 


0.332, 0.498 


0.079 


±0.024 


+0.015 
-0.021 


+0.020 
-0.003 


1.178 ± 0.140 


0.498, 0.664 


0.064 


±0.035 


+0.050 
-0.028 


+0.001 
-0.013 


1.374 ± 0.215 


0.664, 0.830 


0.148 


±0.039 


+0.014 
-0.038 


+0.001 
-0.014 


1.608 ± 0.248 






™obs 

"S 


> 0.75 






-0.830, -0.664 


0.557 


±0.066 


+0.069 
-0.054 


+0.027 
-0.017 


0.758 ± 0.014 


-0.664, -0.498 


0.371 


±0.048 


+0.024 
-0.021 


+0.018 
-0.016 


0.842 ± 0.041 


-0.498, -0.332 


0.258 


±0.046 


+0.034 
-0.028 


+0.017 
-0.023 


0.880 ± 0.053 


-0.332, 0.000 


0.183 


±0.024 


+0.022 
-0.024 


+0.009 
-0.017 


0.914 ± 0.048 


0.000, 0.332 


0.198 


±0.024 


+0.018 
-0.018 


+0.013 
-0.017 


0.922 ± 0.062 


0.332, 0.498 


0.212 


±0.035 


+0.029 
-0.008 


+0.013 
-0.008 


0.892 ± 0.065 


0.498, 0.664 


0.313 


±0.053 


+0.024 
-0.043 


+0.001 
-0.023 


0.885 ± 0.076 


0.664, 0.830 


0.308 


±0.068 


+0.031 
-0.023 


+0.014 
-0.019 


0.831 ± 0.071 



Table 4: Measured cross sections as a function of cos 9* for x° < 0.75 and 
x obs > Q.75. For further details, see the caption to Tabled 
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Figure 1: LO QCD charm-production diagrams, a) direct photon: 7(7 — > cc; 
b) resolved photon: gg — > cc; c) resolved-photon charm excitation: eg — * gc (c in 
proton hemisphere); d) resolved-photon charm excitation: eg — > eg (c in photon 
hemisphere). 
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Figure 2: Differential cross-section da/dx^ for the data (dots) compared with: 
a) various MC simulations (histograms); b) NLO FO predictions after hadroni- 
sation correction (full lines) and at parton level (dashed lines). The inner error 
bars show the statistical uncertainty, while the outer ones show the statistical and 
systematic uncertainties added in quadrature. The jet- energy- scale uncertainty is 
given by the two dashed-dotted lines. In a), each MC distribution is normalised 
to the data, as indicated in the brackets. Also shown in a) is the resolved photon 
distribution (hatched) of PYTHIA and in b) the uncertainty of the NLO prediction 
after hadronisation correction (shaded). In b) the two highest £° bs bins have been 
combined. 
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Figure 3: Differential cross-section da/dx° bs for the data (dots) compared with: 
a) PYTHIA and HERWIG MC simulations (histograms); b) CASCADE (short- 
dashed lines) and NLO FO predictions after hadronisation correction (full lines) 
and at parton level (long-dashed lines). The inner error bars show the statistical 
uncertainty, while the outer ones show the statistical and systematic uncertainties 
added in quadrature. The jet-energy-scale uncertainty is given by the two dashed- 
dotted lines. In a), each MC distribution is normalised to the data, as indicated 
in the brackets. Also shown in a) is the resolved photon distribution (hatched) 
of PYTHIA and in b) the uncertainty of the NLO prediction after hadronisation 
correction (shaded). 
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Figure 4: Differential cross-sections da/d | cos 9*\ (dots) compared with: a-b) 
PYTHIA and HERWIG MC simulations (histograms); c-d) CASCADE (short- 
dashed lines) and NLO FO predictions after hadronisation correction (full lines) 
and at parton level (long-dashed lines). Results are given separately in a,c) for 
samples enriched in resolved photon events and in b,d) for samples enriched in 
direct photon events. The inner error bars show the statistical uncertainty, while 
the outer ones show the statistical and systematic uncertainties added in quadra- 
ture. The jet- energy -scale uncertainty is given by the two dashed-dotted lines. In 
a-b), each MC distribution is normalised to the data, as indicated in the brackets. 
Also shown in c-d) are the uncertainties of the NLO prediction after hadronisation 
correction (shaded). 
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Figure 5: Differential cross-sections da/dcos8* (dots) compared with: a-b) 
PYTHIA and HERWIG MC simulations (histograms); c-d) CASCADE (short- 
dashed lines) and NLO FO predictions after hadronisation correction (full lines) 
and at parton level (long-dashed lines). Results are given separately in a,c) for 
samples enriched in resolved photon events and in b,d) for samples enriched in di- 
rect photon events. The inner error bars show the statistical uncertainty, while the 
outer ones show the statistical and systematic uncertainties added in quadrature. 
The jet-energy-scale uncertainty is given by the two dashed-dotted lines. In a-b), 
each MC distribution is normalised to the data, as indicated in the brackets. Also 
shown as shaded areas in a) and b) are the contribution of the direct photon process 
in PYTHIA to the resolved- enriched sample and the contribution of the resolved 
photon process to the direct- enriched sample, respectively. The uncertainties of the 
NLO prediction after the hadronisation correction are shown as the shaded areas 
in c) and d). <->-, 



